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13  AB&TRACT- 


An  analytical  model  for  the  dynamic  response  of  a  bistable  fluid  amplifier  with 
straight  walls  is  presented.  The  model  represents  the  dynamic  changes  which  take 
place  in  an  amplifier  during  the  early  phases  of  switching.  The  state  of  the  ampli¬ 
fier  at  any  instant  of  time  is  identified  by  a  set  of  system  equations  whose  input 
function  is  the  mass  flow  rate  from  the  control  port.  The  system  variables,  constant! 
and  equations  are  arranged  in  nondimens ional  form  in  order  to  organize  and  facilitate 
the  solution. 

In  the  second  half  of  the  report,  the  system  equations  are  used  to  prepare  an 
analog  simulation  of  the  amplifier.  The  coupling  between  system  variables  and  the 
feedback  mechanisms  present  in  the  model  can  be  easily  identified  from  this  simulatioi 

Computer  results  were  obtained  for  this  simulation  for  a  number  of  input  con¬ 
ditions  and  amplifier  geometries.  One  of  the  more  interesting  results  of  this  study 
is  the  dynamic  change  of  pressure  within  the  separation  region  (bubble).  With  a 
sudden  increase  in  control  flow,  the  bubble  undergoes  an  initial  decrease  in  pres¬ 
sure  followed  by  its  expected  rise  to  a  new  equilibrium  condition.  Because  of  this 
result,  the  jet  attachment  point  moves  toward  the  power  nozzle  for  a  short  duration 
of  time. 
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ABSTRACT 


An  analytical  model  for  the  dynamic  response  of  a  bistable  fluid 
amplifier  with  straight  walls  is  presented.  The  model  represents  the 
dynamic  changes  which  take  place  in  an  amplifier  during  the  early 
phases  of  switching.  The  state  of  the  amplifier  at  any  instant  of 
time  is  identified  by  a  set  of  system  equations  whose  input  function 
is  the  mass  flow  rate  from  the  control  port.  The  system  variables, 
constants  and  equations  are  arranged  in  nond imens ional  form  in  order 
to  organize  and  facilitate  the  solution. 

In  the  second  half  of  the  report,  the  system  equations  are  used 
to  prepare  an  analog  simulation  of  the  amplifier.  The  coupling  between 
system  variables  and  the  feedback  mechanisms  present  in  the  model  can 
be  easily  identified  from  this  simulation. 

Computer  results  were  obtained  for  this  simulation  for  a  number 
of  input  conditions  and  amplifier  geometries.  One  of  the  more  inter¬ 
esting  results  of  this  study  is  the  dynamic  change  of  pressure  within 
the  separation  region  (bubble).  With  a  sudden  increase  in  control  flow, 
the  bubble  undergoes  an  initial  decrease  in  pressure  followed  by  its 
expected  rise  to  a  new  equilibrium  condition.  Because  of  this  result, 
the  jet  attachment  point  moves  toward  the  power  nozzle  for  a  short  dur¬ 
ation  of  time. 
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I. 


INTRODUCTION 


Many  control  systems  which  are  in  use  today  depend  on  an  assort¬ 
ment  of  pneumatic  and  hydraulic  devices.  These  components  utilize  the 
characteristics  of  flowing  fluids  for  actuation,  decision  making  and 
amplification.  The  major  class  of  such  devices  have  moving  parts  and 
depend  on  the  interaction  between  mechanical  and  fluidic  systems.  Some 
common  amplifiers  with  moving  parts  are; 

1.  spool  valves 

2 .  flapper  valves 

3.  ball  valves 

4.  diaphragm  valves 

During  the  past  tweive  years  a  second  class  of  fluid  amplifiers 
has  been  developed.  These  devices  have  no  moving  parts  and  are  referred 
to  as  pure  fluid  amplifiers.  Their  input  is  a  fluidic  signal  whicn  pro¬ 
duces  an  amplified  change  in  the  output  fluid  flow.  The  geometry  of 
the  amplifier  remains  fixed  during  its  operation.  These  pure  fluid  amp¬ 
lifiers  can  be  further  subdivided  into  two  categories.  Those  having  a 
continuous  input-output  relationship  are  called  analog  or  proportional 
amplifiers.  The  second  group  have  digital  characteristics  and  are  known 
as  bistable  or  two-state  amplifiers. 

This  second  class  of  pure  fluid  amplifiers  have  one  or  more  suable 
output  states  and  are  used  in  the  design  of  logic  systems.  Several  of 
these  fluid  logic  devices  are  connected  to  achieve  different  logic  or 
switching  functions.  Four  of  the  more  common  elements  are: 

1.  Induction  amplifier 

2.  Turbulence  amplifier 

3.  Focused  jet  amplifier 

4.  Wall-attachment  amplifier 

Modeling  of  the  fourth  element,  the  wall-attachment  amplifier,  is  the 
subject  of  this  report. 

The  development  of  the  wall-attachment  amplifier  was  first  an¬ 
nounced  at  Diamond  Ordinance  Fuze  Laboratories  (now  Harry  Diamond  Labora¬ 
tories,  HDL)  in  1960.  Operation  of  this  amplifier  depends  on  the  Coanda 
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effect  which  is  the  natural  attachment  of  a  ftee  jet  to  an  adjacent 
wall.  A  low  pressure  region  or  bubble  exists  between  the  wall,  the 
jet  and  the  point  of  attachment.  The  pressure  difference  across  the 
jet  establishes  the  curvature  of  the  jet  as  it  exits  from  tht  output 
duct  on  the  attached  side.  The  addition  of  control  flow  into  the  low 
pressure  region  causes  the  bubble  to  grow.  This  growth  results  in  the 
jet  attaching  to  the  opposite  wall  and  issuing  from  the  output  duct  on 
that  side.  This  process  takes  place  in  a  few  milliseconds  and  estab¬ 
lishes  the  bistable  relationship  for  the  amplifier. 

Initial  studies  on  the  reattachment  of  free  jets  to  inclined  and 
offset  walla  were  done  by  many  different  investigators. i  The  work  of 
Bourque  and  Newman  [1]  has  served  as  the  basis  for  n.uch  of  the  recent 
studies  [2,  3,  4).  The  application  of  this  work  to  the  amplifier  re¬ 
sulted  in  studies  on  steady  state  performance  and  the  effects  of  amp¬ 
lifier  geometry  on  attachment.  The  works  of  Sher  [5]  and  Levin  and 
Manion  [6]  are  indicative  of  these  types  of  studies  [7,  8,  9], 

The  next  phase  of  research  dealt  •  the  switching  process  and 
its  characteristics.  A  number  of  paperi  ave  been  written  or.  the 
subject  [10,  11,  12,  13,  14],  Kirshner  [15]  defined  the  three  types 
of  switching  (1)  terminated  wall  ' (2)  contacting-both-walls  and  (3) 
splitter  switching,  while  Muller  [11]  described  a  critical  attachment 
angle  type  of  switching. 

Most  recent  studies  have  concentrated  on  the  dynamic  modeling  of 
the  switching  process  itself.  Savkar,  Hansen  and  Keller  [14]  reported 
some  experimental  values  for  switching  times  while  Lush  [16,  17]  com¬ 
pared  experimental  values  of  switching  times  with  a  theoretical  pre¬ 
diction.  Recently  theoretical  and  experimental  results  for  the 
contacting-both-walls  type  of  switching  have  been  obtained  by  Ozg^i 
and  Stenning  [22].  Although  many  investigators  [18,  19,  20,  21]  have 
contributed  to  the  study,  a  complete  description  of  the  switching  pro¬ 
cess  has  not  been  obtained.  This  project  was  undertaken  in  an  attempt 
to  provide  some  understanding  of  the  early  phases  of  dynamic  switch¬ 
ing  in  the  wall-attachment  amplifier,  , 


II.  MODELING  THE  DYNAMIC  RESPONSE 


A.  Amplifier  Description 

I 

The  bistable  amplifier  to  be  analyzed  is  the  standard,  sloping, 
straight  wall  type.  The  control  port  is  perpendicular  to  the  power 
jet  inlet  and  has  width  C.  The  geometry  of  this  particular  ampli¬ 
fier  is  shown  in  Figure  1.  The  amplifier  is  also  considered  to  be 
symmetric  about  its  centerline  with  the  parameters  D,  a  and  de¬ 
fining  the  offset,  wall  angle  and  distance  to  the  vents. 

In  an  equilibrium  state  the  jet  is  attached  to  one  wall  of  the 
amplifier  in  such  a  manner  as  to  form  a  region  of  uniform  pressure 
called  the  separation  bubble.  This  bubble  region  maintains  a  con¬ 
stant  volume  and  mass  as  long  as  the  net  mass  flow  into  the  region  is 
zero.  The  centerline  of  the  jet  is  as6umed  to  be  a  circular  arc  of 
radius  R  which  leaves  the  center  of  the  power  jet  at  an  angle  g  . 

This  arc  intersects  the  amplifier  wall  at  an  angle  0  and  at  a  point 
determined  by  the  geometry. 

Outer  boundary  of  the  separation  bubble  is  formed  by  the  attaching 
streamline  of  the  fluid  jet  profile.  This  arc  is  given  by  the  distance 
y  from  the  circular  arc  measured  along  a  radius  vector.  The  remaining 
boundary  of  the  bubble  is  formed  by  the  fixed  boundary  of  the  amplifier 
The  jet  attaches  to  the  wall  at  a  point  defined  by  Xf.  Geometry  of  the 
amplifier  along  with  the  attached  jet  is  shown  in  Figure  2. 

Before  an  analytical  model  can  be  developed  several  assumptions, 
must  be  made  concerning  the  jet  behavior  and  the  other  interactions 
within  the  amplifier. 

1.  The  attached  jet  has  a  centerline  which  follows  a  circular 
arc  of  radius  R. 

2.  The  attached  jet  is  turbulent  and  has  a  Goertler  velocity 
profile. 

3.  The  thickness  of  the  two  dimensional  jet  is  d. 

4.  Pressure  inside  the  separation  bubble  is  uniform  throughout 
the  region. 
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5.  The  width  of  the  jet  is  small  compared  to  the  radius. 

6.  Momentum  interaction  between  the  control  and  power  jet  occurs 
at  a  discrete  point. 

7.  The  fluid  within  the  separation  region  behaves  as  an  ideal  gas. 

From  these  assumptions  and  the  geometry  of  the  amplifier  a  dynamic 
model  of  the  amplifier  can  be  developed.  This  model  which  represents 
the  stare  of  the  amplifier  at  any  instant  of  time  will  be  represented 
by  a  set  of  system  variables.  These  variables  are  related  by  a  set  of 
algebraic  and  differential  equations  referred  to  as  system  equations. 

B.  Ceometric  Relations 

A  number  of  system  equations  can  be  determined  from  the  geometry 
of  the  amplifier  shown  in  Figure  2.  By  equating  distances  perpendic¬ 
ular  to  the  wall  of  the  amplifier  the  following  equation  results. 

R  cos  (  a  +  3  )  =  r  cos  9  +  (D  +  ^  -  C  tan  a  )  cos  a  (1) 


At  the  attachment  point  the  length  r  can  be  written, 

X 

r  =  R  -(£)  R(  S  +  a  +  0  )  -  n(^  ) 


(2) 


where  n  defines  the  attaching  streamline,  o  is  the  Goertler  spread 
parameter  and  Xq  determines  the  virtual  origin  of  the  Jet.  The  small 
angle  $  can  be  determined  from  the  other  dimensions. 


8  in  4> 


(b/2)  sin  B _ 

?  i  /? 

(i  -  bR  cos  9  +  b  /4) 


(3) 
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By  introducing  a  new  variable  T  } 


0  +  8  +  a  +  9 


(4) 


the  expression  for  the  bubble  volume  (V  )  can  be  simplified.  The  de- 

D 

velopment  of  the  following  equation  is  given  in  Appendix  A. 


d  {  \  (R  -  |  cos  (  6  +  0  ))2  r  -  (R  -  |  cos  (  8  +  0  ))J  R \ 


R  .  2  T  1  2 


+  (  H  —  )  g  -  j  r  cos  ®  [tan  (  8  +  0  +  ot  )  +  tan  9 

1 


+  ilD.ctana!2  f  B-tifeJL  ,<*  1 

+  2  ID  c  tan  a  j  l  sin(9Q  _  Q  g  >  $  ) 


(5) 


Finally,  the  point  of  reattachment,  Xr>  can  be  calculated. 


*»  [R  sin  8  -C]/cos  a  +  (r  cos  9  )  tan  a  +  r  sin  0 


(6) 


C.  Fluid  Relations 

Additional  equations  can  be  obtained  by  considering  the  conditions 

of  the  fluid  within  the  bubble  region.  Assuming  that  the  entrapped 

fluid  behaves  as  an  ideal  gas  with  uniform  pressure  Pn,  volume  VB  and 

mass  M  the  following 
B 

PBVB  "  W 


(7) 


represents  the  equation  of  state  where  R^  (gas  constant)  and  T  (temp¬ 
erature)  remain  constant.  Mass  is  being  added  to  the  region  from  the 
control  port  at  the  rate  mc>  The  power  jet  entrains  mass  from  the  re¬ 
gion  at  a  rate  m  and  recirculates  mass  at  a  rate  m  .  The  change  of 
e  r 

mass  inside  the  bubble  can  be  expressed  as. 


dt 


m  +  m  -  m 
ere 


(8) 


The  angle  6  can  be  determined  by  considering  the  momentum  interaction 
between  the  control  and  power  jets, 


M  C 


where  M  is  the  mass  flow  rate  of  the  power  jet. 

The  motion  of  the  power  jet  is  represented  by  a  term  which  de¬ 
pends  on  the  rate  of  change  of  radius  R  in  the  following  curvature 
equation. 


3M  Ao  dR  _  M2/  c  bd2  _ 

bda  Y  Rf  dt  R  B 


(10) 


Development  of  this  equation  along  with  a  discussion  of  the  additional 
term  is  given  in  Appendix  B.  Parameter  P  is  the  density  of  the  power 
jet  and  is  assumed  to  remain  constant. 

Using  the  Goertler  profile  for  the  velocity  distribution  of  the 
jet,  the  entrainment  and  recirculation  flow  rates  can  be  determined. 

See  Appendix  C. 
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The  model,  at  this  point,  describes  the  dynamic  growth(collapsing) 

of  the  separation  region  (bubble)  due  to  changes  in  control  mass  flow 

m  .  With  m  as  the  input,  the  state  of  the  amplifier  at  any  instant 
c  c 

of  time  is  described  by  thirteen  variables  listed  in  the  table  below. 


Table  1 


Notation_ Description 


. . —  — . -  *  ■-  -  —  -  "  — 1 

1.  m 

entrained  mass  flow  rate 

2.  m 

recirculated  mass  flow  rate 

3.  m 

control  mass  flow  rate  (input) 

4. 

bubble  mass 

3- 

bubble  pressure 

6.  Xr 

reattachment  length 

7.  R 

radius  of  curvature  of  jet 

8.  r 

radius  to  reattachment  point 

9‘  VB 

bubble  volume 

10.  3 

initial  jet  deflection  angle 

11.  0 

angle  of  attachment 

12.  4) 

small  angle 

n.  r 

' 

summation  of  four  angles 

With  the  control  mass  flow  rate  specified,  the  system  equations  con¬ 
sist  of  the  twelve  equations  numbered  1  through  12. 

A  simplified  version  of  a  nonlinear  block  diagram  for  the  proposed 
model  is  shown  in  Figure  3.  Each  block  represents  a  system  equation 
and  is  numbered  accordingly.  The  output  of  each  block  can  be  deter¬ 
mined  when  each  of  the  indicated  inputs  are  known.  Some  of  the  con¬ 
necting  lines  have  been  omitted  for  simplification. 

The  primary  feedback  paths  involve  m  and  m^.  Since  the  entrained 
mass  flow  rate  is  subtracted  it  has  a  stabilizing  effect  on  the  output 
radius,  R.  Instability  is  caused  by  the  positive  feedback  of  the  re¬ 
circulation  rate,  m^.  See  Equation  8.  Another  feedback  path  results 
due  to  the  effect  of  the  radius  on  the  volume  (block  5)  which  in  turn 
effects  the  bubble  pressure  (block  7).  A  predominant  feedforward  path 

exists  through  m  to  B  to  V  to  P  . 

c  B  B 

The  input  for  the  system  is  the  control  flow  rate  mc  and  repre¬ 
sents  a  known  function  of  time. 


III.  NONDIMENS IONAL  MODEL 

It  is  desirable  to  have  the  modeling  process  be  relatively  inde¬ 
pendent  of  the  exact  dimensions  of  the  amplifier  and  the  fluid  proper¬ 
ties.  For  this  reason,  a  set  of  nondimens ional  system  variables  is 
established  based  on  the  physical  variables.  These  new  quantities  are: 

[R/b],  (r/b) ,  (Vg/b2d],  (Mb/  pb2d],  [mc/M] 

[m./M],  (me/M],  (Pg/^/p  C2d2],  [Xr/b] 

[  e  J>  l  e  l  $  I*  l  r  J 

The  independent  variable,  time,  is  also  made  nondimens ional  by  the 
following  substitution. 
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Substitution  of  these  variables  into  the  system  equations  yield  a  new 
set  of  nondimens ional  equations  with  new  coefficients.  These  coeffi¬ 
cients  can  be  collected  into  a  set  of  constants  which  become  the  non- 
dimensional  parameters  for  the  model.  The  nondimens ional  system  vari¬ 
ables  X^  and  parameters  are  defined  in  Tables  2  and  3. 


Table  2. 

Identification  of 

the  nondimens ional 

system  variables. 

Physical  Variable 

Units 

System  Variable 

1. 

(n»e/M) 

1 

X1 

2. 

(mr/M) 

1 

X2 

3. 

(mc/M) 

1 

(input; 

4. 

CM./  P  b2d) 

B  M2 

1 

X4 

5. 

(V  2,2  } 

1 

X5 

p  C  d 

6. 

<xr/b) 

1 

X6 

7. 

(R/b) 

1 

X7 

8. 

(r/b) 

1 

X8 

9. 

<Vfi/b2d) 

1 

X9 

10. 

6 

radians 

X10 

11. 

e 

radians 

X11 

12. 

$ 

radians 

X12 

13. 

r 

radians 

X13 
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Table  3.  Identification  of  the  system  parameters 


Physical  Constants 

System  Parameters 

2  2  2 

,  RaT  p  C  d  gc 

M2 

K1 

2.  (C/b) 

K2 

3 .  1  -  tanh  ( n  ) 

K3 

4.  (pb2d/M)  «  t/T 

K4 

5.  (3/4  o  ) 

K5 

6.  D/b 

K6 

7.  a 

K7 

8.  (  n  /  a  ) 

K8 

9.  (XQ/b)  -  a  / 3 

K9 

10.  ?J(V2/p  C2d2) 

K10 
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The  nondimens ional  system  equations  in  terms  of  and  are  as  follows: 


X7cos(K7  +  X1Q)  -  XgcosX11  +  (K6  +  j)  cosK?  -  K2sinK? 


(13) 


Xg  -  X?  -  KgX7(X10  +  X11  +  K7)  -  KgKg 
sinX12  “  ({  sinX10)/(X7  -  X7cosX10  +  0.25)l/2 


(14) 

(15) 


13 


xio  +  xn  +  xi2  + 


«  (~)  (X?  -  0.5cos(X10  +  X12)2]  Xn 

X2  X3 

•  1*7  -  7  “s  <X10  +  X12)l  KgX,-^  +  -|1 


-  |  Xg  cos2X11[tan(X13-  Xu)  +  tan  XU1 


(16) 

07) 


1  2 

+  y  (Kg  -  K2  tan  K7)  (ain(XlQ  +  X12)  cos  K? ) 

COS  (K7  +  X10  +  X12) 


Xg  «  [X7  sin  X10  -  K2  ] /cos  K -j  +  Xg  (cos  X^  tan  K7  +  sin  X^]  (18) 


x5  *9  »  *4  h 


(19) 


dx4 
d  T 


x3  +  x2  -  x3 


tan  X1Q  =  X3/K2 


(20) 

(21) 


(k2) 


1/2  dX7 


(K9) 


T72  (X13X7}  ’  TT 


(K2) 

‘  xf  "  X5  '  K10 


(22) 
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Xt  -  (K5X7(X13  -  X12)  +i]  -  ~  (23) 

!  1/2 

X2  -  [K5X7(X13  -  X12)  K3  (24) 

There  are  13  system  variables  which  define  the  state  of  the  ampli¬ 
fier  at  any  instant  of  time.  The  first  three  (X^,  X2,  X3)  represent 
flow  rates  with  X3  serving  as  the  input  to  the  system.  The  bubble  mass, 
pressure  and  volume  are  represented  by  variables  X^,  X^  and  X^.  Three 
distances  within  the  amplifier  are  associated  with  X,,  X,  and  XD.  Fi- 
nally,  the  last  four  system  variables  replace  the  angles  which  describe 
the  bubble  geometry. 

Of  the  10  system  parameters  listed  in  Table  3,  the  only  one  which 
has  units  associated  with  it  is  K^.  Since  this  parameter  relates  the 
model  time  with  physical  time  it  must  have  the  units  of  time.  The 
other  nine  parameters  must  have  specific  values  in  order  to  generate 
a  solution,  however,  is  only  used  to  convert  the  response  time  to 
a  particular  amplifier  and  fluid. 

The  12  system  equations  consist  of  10  algebraic  equations  and 
two  differential  equations.  Most  of  the  equations  are  nonlinear  and 
contain  cross-coupling  between  variables.  For  example,  equations  13,  14 
and  16  contain  Xg,  and  X^g  and  a  direct  solution  for  any  one  of 
the  variables  is  extremely  difficult. 

IV .  ANALOG  MODEL  AND  SIMULATION 

The  analytical  model  for  the  amplifier  can  now  be  translated  into 
the  form  of  an  analog  diagram.  This  diagram,  shown  in  Figure  4,  utilizes 
standard  analog  symbols  which  are  explained  in  Appendix  D.  The  flow  of 
information  and  feedback  paths  which  were  present  in  the  original  block 
diagram  are  also  present  in  the  analog  diagram. 

Each  of  the  system  variables  represents  a  possible  output  and  ap¬ 
pears  somewhere  on  the  diagram.  The  necessary  inputs  to  the  model  are 
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the  values  for  the  system  parameters  the  functional  relation 
for  the  input  X^(  T)  and  the  initial  values  for  variables  X^  and  X^. 

As  a  result,  both  a  steady  state  and  dynamic  simulation  are  necessary 
to  predict  the  response. 

A.  Steady  State  Response  Simulation 

The  steady  state  or  equilibrium  state  of  the  amplifier  is  reached 
when  the  sum  of  all  mass  flows  into  the  separation  bubble  is  zero.  In 
this  case  the  bubble  remains  stationary  and  all  derivatives  of  the 
system  variables  become  zero.  Thus,  the  system  equations  (Equations  13 
through  24)  reduce  to  a  set  of  simultaneous,  nonlinear  algebraic  equa¬ 
tions.  The  steady  state  simulation  consists  of  an  iterating,  FORTRAN 
program  which  solves  these  algebraic  equations  and  predicts  the  equi¬ 
librium  values  for  the  12  system  variables.  The  inputs  for  the  pro¬ 
gram  are  the  control  flow,  Xy  and  the  10  constants  which  depend  on  the 
amplifier  geometry,  fluid  properties  and  approximate  jet  flow  solutions. 

Computer  results  were  obtained  for  values  of  ranging  from  0  to 
0.5.  This  corresponds  to  a  control  mass  flow  rate  from  0  to  507o  of 
the  power  jet  mass  flow  rate. 

In  addition,  two  parameters  and  which  are  the  ratios  of  con¬ 
trol  port  width  and  offset  distance  to  power  nozzle  width  were  assigned 
several  different  values.  Parameters  Ky  Kg  and  Kg  were  also  changed 
to  correspond  to  several  different  values  of  spread  parameter  a  . 

The  results  for  two  cases  are  shown  in  Figures  5  through  14.  The 
steady  state  values  for  each  of  the  system  variables  are  indicated  for 
a  range  of  values  for  control  input  Xy 

B.  Dynamic  Response  Simulation 

For  this  portion  of  the  response  an  analog  simulation  was  sought. 
The  system  equations  and  block  diagram  lend  themselves  to  a  direct 
analog  simulation  with  the  nondimens ional  time  variable  T  functioning 
as  the  single,  independent  variable.  In  addition,  the  dynamic  response 
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of  any  one  of  the  system  variables  might  be  of  interest  and  could  be 
easily  recorded.  However,  because  of  the  large  number  of  nonlinear¬ 
ities  present  in  the  model  an  analog  computer,  capable  of  handling 
the  problem,  was  not  available.  For  this  reason,  an  alternate  means 
of  solution  was  selected, 

A  packaged  subroutine  which  simulates  an  analog  computer  has  been 
developed  by  Lehigh  University  for  use  on  the  CDC  6400  digital  computer. 
The  program,  referred  to  as  LEANS,  was  used  to  provide  the  analog  sim¬ 
ulation  for  this  problem.  It  employs  the  block  diagram  technique, 
commonly  considered  the  first  step  of  any  analog  simulation.  The  pro¬ 
gram  contains  a  full  array  of  integrators,  summers,  constants,  mult¬ 
ipliers  and  other  nonlinear  elements  which  can  be  coupled  in  any  man¬ 
ner.  Since  the  digital  computer  uses  floating  point  arithmetic,  time 
and  magnitude  scaling  are  unnecessary;  and  the  us jr  may  program  di¬ 
rectly  from  the  block  diagram. 

The  LEANS  simulation  was  prepared  directly  from  Figure  4.  The 
system  equations  were  used  to  provide  the  particular  functional  re¬ 
lationships  that  exist  between  the  system  variables.  The  required 
inputs  to  this  analog  simulation  are  the  numerical  values  for  the  10 
constants,  initial  conditions  for  X^,  X^  and  X^  and  the  final  value 
for  X^ .  The  input  was  programmed  to  change  exponentially  from  its 
initial  to  final  value.  This  type  of  input  would  commonly  occur  in 
practice  and  has  a  much  faster  rise  time  than  the  dynamic  growth  of 
the  bubble. 

Results  of  this  portion  cf  the  study  are  shown  in  Figures  15 
through  36,  The  initial  value  for  each  system  variable  (  T  =0)  and 
the  final  value  (  T  »  00  )  must  correspond  to  those  values  predicted 
by  the  steady  state  simulation.  For  example,  in  Figure  20  the  initial 
value  for  is  5.16  and  the  final  value  is  6.3  as  changes  from 
0.10  to  0.14.  Compare  these  values  to  those  shown  in  Figure  9  for 
of  0.1  and  0.14. 
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V. 


RESULTS 


With  the  analog  simulation  on  the  computer  it  became  possible 
to  make  a  series  of  parametric  studies.  Changes  in  amplifier  geom¬ 
etry,  fluid  conditions  and  input  specifications  could  be  introduced 
by  merely  changing  the  parameters  of  the  simulation.  Steady  state 
results  are  shown  in  Figures  5  to  14  for  two  different  cases.  Figures 
15  to  36  contain  the  results  of  the  dynamic  simulation  for  two  cases 
A  and  B.  The  parameter  values  for  cases  A  and  B  are  given  in  Figures 
15  and  26  respectively. 

A.  Steady  State  Results 

Entrainment  and  recirculated  mass  flow  rates  (X^,  X^)  appear  in 
Figure  5  to  be  linearly  dependent  on  the  control  mass  flow  input,  X^> 
but  independent  of  parameter  .  Bubble  pressure,  X,.  in  Figure  7, 

Jet  radius  of  curvature,  X^  in  Figure  9,  and  attachment  radius,  X^ 
in  Figure  10,  level  off  as  the  input  approaches  0.5.  For  the  given 
range  of  control  flow,  the  reattachment  angle  6  decreases  at  first 
and  then  begins  to  increase.  With  low  control  flow  the  bubble  is 
very  small  and  the  jet  reattaches  close  to  the  supply  nozzle.  As  a 
result,  the  reattachment  angle  must  be  somewhat  larger. 

B.  Dynamic  Results 

Results  were  obtained  from  the  computer  simulation  for  both  large 
and  small  changes  in  input  signal  X^,  the  ratio  of  mass  flow  from  the 
control  port  to  supply  flow.  The  exponential  input  signal,  shown  in 
Figures  15  and  26,  reaches  its  final  value  very  rapidly  (  T  "3). 

The  dynamic  changes  which  take  place  in  the  system  variables 
for  the  small  signal  input  are  shown  in  Figures  16  to  25.  The  most 
significant  result  is  the  change  in  bubble  pressure  (X^)  shown  in 
Figure  18.  The  model  predicts  a  sudden  decrease  in  pressure  as  the 
control  signal  is  applied.  The  bubble  pressure  recovers,  overshoots 
and  finally  settles  down  to  its  new  equilibrium  value  at  T  =  15. 

The  reason  for  this  inverse  response  (output  variable  initially  de¬ 
creasing  with  an  increase  in  input  signal)  is  the  feedforward  loop 


-15- 


Vi'.  iW  ii’.iM'i  diagram  w'-ficl  (J :  gur*.  3).  An  ii>l  f«aj«  >n  Ir.f  input  sig¬ 
nal  feeds  throvg'-  cl or**  •»,  1,  4  an;l  3  an«l  Has  c  .decreasing  effect  on 
ike  output  variables.  The  physical  fcXj-lai.a'.  ion  it  thar ,  although  an 
sde  it  ion  of  control  flow  increases  the  mass  in  toe  bubble  which  in- 
» teases  the  bubble  pressure,  it  has  a  more  direct  effect  of  increas¬ 
ing  the  bubble  volume.  Duo  to  the  compressibility  ol  air  Inside  the 
hubi.lt  a  decrease  in  bubble  pro»sure  remits.  This  invirce  response 
iii'.  -  appears  in  toe  radius  oi  airvatmv  lor  in*-  jet  (Figure  20)  and 
ti..'  rest  v  achmei't  radius  in  F'pure  21.  In  both  cases,  however,  there 
it  no  oversJ'.i.ot  i ng  of  the  1  irt.l  value. 

The-  ci'.t-i  a:,  this  res  pens*.  ciu’.o  he  obuerci-J  lhro<j|<»  tire  r«  at  - 
t  .i.'.h'jer.c  point  >hr  i-r.  in  ?  i  sure  it.  Y'\  a  snort  cerieV.  o*  tine,  the 
teat  Cat  hment  uo  i  ne  move*,  very  slightly,  toward  the  supp:;-  nozzle. 

T.-v  hu/'ls  .  It- me  in  figure  22.  on  the  other  hand,  increases  continu¬ 
ously  from  its  .nit.ai  equilibrium  value-  to  its  final  value. 

fer  this  narticula.  ch..'*nitc  m  control  tluw  the  reat taefment  angle  i* 
is  vii  'veil  to  a  lover  equilibrium  value  (see  Figure  13).  The  ar-sle, 
however,  increases  1  jt  a  rbort  period  of  time  before*  jtcreasmg  to  its 
final  value.  The  dynamic  rer.por.se  of  this  variable,  35  veil  a»  the 
others,  take?  abo.it  15  units  of  "  to  reach  its  r.«v  equilibrium  position, 
it  the  value  of  K/  (see  Tarlc  3).  is  known  for  a  particular  amplifier, 
then  the  tccal  time  for  f  lit  response  becc-cws  15  K.  in  the  t  ti.ir  units  ci 
K.  . 

The  ayiianiic  respcr.E  t-  of  the  .t.'ulation  to  a  large  cliange  in  con¬ 
trol  flow  vas  also  sougnt  .  The-  sudden  change  in  control  flow  variable 
Y.  from  0.1  to  0.3  lad  the  form  sh-run  tit  Figure  ,10.  Due  tc  the  ab¬ 
sence-  wi  iide  vents  and  t  splitter  the  jet  rema  is  attached  .o  the 
same  wall  but  enclosing  a  much  larger  separation  rtgiop.  T:<e  final 
equilibrium  position  for  the  amplifier  is  reached  at  about  T  ■  30. 

The  equilibrium  values.  for  each  i  i  tnc  system,  variables  louche*,  those 
erven  by  the  steady  stall  £  »:<-u>al  l  c-u  i  it  >1.  -  0 .  i  ai.c  the  cor  responding 
pau^'tet  values  tor  *-•'  . 
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The  dynamic  responses  of  each  of  the  system  variables  for  this 
case  are  given  in  Figures  27  through  36.  The  same  trends  are  evi- 

I 

dent  for  this  case  as  the  small  input  case  but  seem  more  pronounced. 
Bubble  pressure,  for  example,  undergoes  a  much  larger  decrease  ini¬ 
tially  as  shown  in  Figure  29.  The  inverse  response  is  still  present 
in  variables  X^,  Xy  and  Xg  as  displayed  in  Figures  30,  31  and  32. 

For  the  large  value  of  control  flow  (Xg  =  0.3),  the  final  value 
for  attachment  angle  is  slightly  larger  than  the  initial  value  (see 
Figure  13).  However,  during  the  dynamic  response  the  angle  overshoots 
the  value  rather  significantly  as  displayed  in  Figure  35.  Once  again 
it  takes  30  units  of  T  for  the  variable  to  come  to  equilibrium. 

VI.  CONCLUSIONS 

A  nonlinear,  dynamic  model  for  the  bubble  growth  within  a  wall 
attachment  amplifier  has  been  presented.  Using  nond imens ional  system 
variables,  this  model  was  described  in  analog  form.  From  a  computer 
simulation,  the  response  to  a  sudden  change  in  control  mass  flow  has 
been  studied.  These  results  produced  the  following  list  of  conclusions. 

1.  The  separation  bubble  undergoes  an  initial  decrease  in  pressure 
followed  by  an  overshoot  of  its  final  value.  If  the  separation 
region  should  encounter  a  vent  during  this  overshooting  period  no 
mass  would  flow  from  the  region  through  the  vent. 

2.  For  a  short  period  of  time  the  reattachment  point  moves  toward  the 
supply  nozzle,  before  moving  outward, 

3.  The  reattachment  angle  undergoes  a  rapid  increase  and  an  overshoot 
1  of  its  final  value.  During  this  overshoot  period  the  amplifier 

cou’d  possibly  submit  to  a  "critical  angle"  type  of  switching  even 
though  the  final  value  for  the  angle  does  not  predict  it. 

I 

4.  The  time  it  takes  for  the  amplifier  to  reach  equilibrium  (response 
time)  is  longer  for  a  large  change  in  control  flow  than  for  a  small 
change. 


The  simulation  results  are  presented  in  nondiraensional  form 
but  with  the  fixed  dimensions  of  the  amplifier  known,  the  vari¬ 
ables  can  be  converted  to  actual  units.  A  value  for  the  par¬ 
ameter  will  convert  the  time  scale  to  appropriate  time  units. 


APPENDIX  A 


The  volume  of  the  bubble,  as  shown  in  the  figure  below,  is  the 
product  of  the  amplifier  thickness  d  and  the  area  enclosed  by  curve 
FGH.  This  area  can  be  calculated  by  computing  the  large  area  bound 
by  the  curve  EGH,  subtracting  the  triangular  portion  EIH  and  adding 
the  small  area  inside  curve  FIG. 


Area  enclosed  by  EGH  is  denoted  by  and  can  be  determined  as 
follows : 

Aj  -  /  /  p  dP  d  A 

where  A  and  p  are  the  polar  coordinates  of  a  point  measured 
from  point  E  and  line  EG. 

Thus , 

A  =  f  i  p2  d  A 
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and  from  the  Coertler  profile  assumption: 

p  -  R-y 

y  *  n(Xj  +  x)  -  yt  + 


Substituting 


0  +  a  +  8+,0  2 

*1  »  /  Cj  )  [R-^R  X  —  y i  1  d  X 

0 


0  +  a+8  +  €^ 


A l  -  f  (p  l(R  -  yx)  -  2(R  -  yx)  (~)  X  +  (r—^)  )d  X 


2  „  3 


Ai  c  4>  (R  -  Xi>  r  -  (R  "  >1>  “a^  -T*  +  (JI5,£)  “T 


where 


0  +  a  +  0  +0 
(-j)b  cos  (3+0  ) 


The  areas  of  the  two  large  triangles  are  and  A^  and  can  be  easily 
calculated. 

■  y  (r  sin  0  )  (r  cos  0  ) 

A3  "  "2  lr  coa  (  ®  )  tan  (0+6  +  cx  )  ]  (r  cos  0  ] 


The  small  triangle  FIG  has  area  A^  where 


A  o  —  Id  -  C  tan  a  ]  \s—  UL  +  ■  t 1  .s>  (9.°  ±  a- 
a4  2  c  tan  a  j  iflln  (90  _  a  .  p  .  0  ) 
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Finally,  the  bubble  volume  can  be  calculated  combining  the 
four  areas. 


VB  ”  d^Al  A2  “  A3  +  A4^ 


-<■(  i,R-yi,2  +(^)2-f-3 

12  2 

■  2  (r  cos  0  (tan  0  +  tan(  $  +  S  +  a  ) )  ] 

+?  '»•««—  i2  i-Ksf9o  - i  iT-Vr}>> 
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APPENDIX  B 


The  relationship  between  power  jet  radius,  R,  and  the  pressure 
difference  across  tne  jet  can  be  obtained  by  writing  the  momentum 
equation  in  the  r  direction. 


a  V  a  V  V .  -  V  V  . 

9  r  .  v  jr  .  _iJ I  _  _J> 

3  t  Vr  Tr^  r  3  4>  r 


Assuming  V  is  independent  of  4>  and  t. 


3  Vt  .  li  =  zJL  3  P 
r  3  r  r  p  3  r 


substituting. 


r  ■  R  +  y 


vr  ■  3?  +  Vy 


3  P 

p  3r 


one  obtains  the  following. 


(15 

vdt 


v  ) 
y 


3Vv 

3  y 


V.2 

_ 


(R  +  y) 


^1  3  P 

p  3  y 


noting  that  ^  and  R  >  >  y  and  neglecting  terms  the  following 

equation  results. 


V  V 

15  (  2  y,  .  x  .  1  — M 

dt  y  '  R  p  3y 
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Integrating  with  respect  to  y  from  +  °°  to  -  ®  , 


( 


dy) 


dR 

dt 


(P*  -  p  ) 

15  00 


where  J  is  the  power  jet  momentum,  P  is  the  bubble  pressure  and  P^, 
s  B 

is  the  pressure  on  the  '■■Dposite  side  of  the  jet.  This  is  the  familiar 

dR 

jet  curvature  equation  with  an  additional  term  that  depends  on  — . 

dR 

For  the  steaoy  state  case  =  0)  it  reduces  to  the  simpler  curva¬ 
ture  expression.  The  integral  can  be  completed  with  this  result. 


+  OO  T 

pv  ]  &  -  i. 

y  dt  R 


p  -  p 
B  * 


Assuming  that  velocity  V 


y 


has 


the  Goertler  distribution, 


3  P  Js 

f  - - 

1  16XJ 


1/2 

} 


JL_D _ 

(en  -  e-n)2 


tanh  n 


+  oo 

} 

-OO 


dR 

dt 


a  differential  equation  is  obtained  which  has  a  coefficient  that  de¬ 
pends  on  the  position  along  the  jet  centerline,  X.  The  supply  jet 
velocity  is  noted  by  Ug . 


s  J  \  dR  3s  r 
o  X  dt  ‘  R  B 


The  average  value  of  this  coefficient  will  be  determined  as  X  varies 
from  zero  to  its  maximum  value  of  R  T  , 

Caxy  = 

1 _ 

r  r 
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Ipi 

M  a 


where  T  is  the  total  angle  swept  out  by  the  bubble.  Using  this 
expression  and  writing  the  supply  Jet  velocity  and  momentum  in  terras 
of  the  mass  flow  from  the  supply  nozzle  (M),  the  following  equation 
is  obtained. 


'v 

—2—  (<£> 
r  r  w 


"2  <'R^ 


-  P 

B  00 
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APPENDIX  C 


The  bubble  region  shown  in  the  figure  below  is  assumed  to 

have  a  uniform  pressure.  Mass  is  being  added  to  the  region  from 

the  control  port  at  a  rate  m^.  The  main  jet  entrains  mass  from 

the  region  at  a  rate  m  and  recirculates  mass  at  a  rate  m  . 

e  r 


m 

c 


(1)  Entrained  mass  flow  rate 


The  total  entrainment  at  some  point  X  is  obtained  by  calculating 
the  total  mass  flow  on  the  bubble  side  of  the  jet  and  subtracting  one 
half  of  the  inlet  flow  rate. 


j.b. 


m  =  p  d  /  udy 
ex  JQ 


0.5  M 

Using  the  Goertler  distribution. 


m 


ex 


00  2 

.  ,  r  ,3M  o 
P  d  j  ( - 5 — 5-) 

0  4  p  Xd  b 


1/2 


where  j.b.  indicues  the 
jet  boundary 


X  2 

- -  sech  n  d  n  -  0.5  M 


m 


ex 


ex 


1/2 


M (-^ ^ — )  /  sech  n  d  n  -  0,5  M 

1/2 
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The  total  entrainment  along  the  entire  edge  of  the  bubble  is  de¬ 
sired,  Substituting, 

R(  8  +  a  +  0  )+x 

o 

g  b 
3 


results  in  the  following: 


M( 


3R 


1/2 


4b  0 


(  8 


+  OL  +  0  )  + 


-  0.5  M 


(2)  Recirculated  mass  flow  rate 


The  total  amount  of  mass  that  recirculates  into  the  bubble  can 
be  calculated  by  integrating  the  velocity  profile  between  the  attach¬ 
ing  streamline  (  n  )  and  4-  00  at  the  attachment  point.  This  point 
is  defined  by  this  attaching  streamline  and  the  boundary  of  the  amp¬ 
lifier. 


P  /  /  udyd2 

Y  „  0  y  2 


p  d  /  (■ 


3M2X 


1/2 


2  ,2V  , 

p  d  b  g  4 


sech  n  d  n 


M( 


3X 


1/2 


4  0  b 


(1  -  tanh  H ) 


Substituting  for  X 


1  /2 

m  =  fal>  ~  ■  (  8  +  a  +  0  )  +  7]  (1-  tanh  n ) 

r  hd  o  4 
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APPENDIX  D 


This  appendix  contains  a  list  of  the  standard  analog  elements 
and  their  corresponding  mathematical  operations. 


Element 


Symbol 


Mathematical 
Operat ion 


Integrator 


t 

2  ■  /  <X+Y)dt 

0 


Multiplier 


Nonl inear 
Funct ion 


X 

m 

Y 


Z  »  XY 


Z  -  f(X,  Y) 
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Figure  1.  Geometry  of  the  Bistable  Amplifier 


Figure  2.  Geometry  of  the  Attached  Jet. 
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Cases  1  and  2 


Figure  10.  Steady  State  Results  for  Reattachment 


Nondimensional  Control  Flow  X 


Figure  16.  Dynami 


Figure  17.  Dynamic  Response  of  Bubble  Hass  (X.)  for  Case 


Figure  18.  Dynamic  Response  of  Bubble 


Nond imensional  Time 


Figure  22.  Dynamic  Response  of  Bubble  Volume  (X^)  for  Case  A. 


Figure  24,  Dynamic  Response  of  Angle  0  (X^)  for  Case  A 
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Figure  26.  Dynamic  Input  for  Case  B  -  Large  Change  in  Control  Flow  (X  ). 
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igure  28.  Dynamic  Re 


Figure  29.  Dynamic  Response  of  Bubble  Pressure  (X^)  fot  Case 


Figure  30.  Dynamic 
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Figure  34.  Dynamic 
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Figure  35.  Dynamic 


